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Measurements of neutrino oscillations using the disappearance of muon neutrinos from the Fer- 
milab NuMI neutrino beam as observed by the two MINOS detectors are reported. New analysis 
methods have been applied to an enlarged data sample from an exposure of 7.25 x 10'^° protons on 
target. A fit to neutrino oscillations yields values of |Am^| = (2.32to og) x 10~^ eV^ for the atmo- 
spheric mass splitting and sin^(2^) > 0.90 (90% C.L.) for the mixing angle. Pure neutrino decay and 
quantum decoherence hypotheses are excluded at 7 and 9 standard deviations, respectively. 
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PACS numbers: 14.60.Pq, 14.60.Lm, 14.60.St, 29.27.-a, 29.40.-n, 29.40.-n, 29.40.Mc 



Neutrino masses and flavor mixing influence the role of 
neutrinos in fundamental physics processes [l[ and may 
point to the mechanism that gives rise to the matter- 
antimatter asymmetry in the Universe 0- A variety of 
phenomena observed with neutrinos originating in the 
Earth's atmosphere or the Sun and those produced by 
nuclear reactors or accelerators can be described consis- 
tently by quantum-mechanical mixing of the weak flavor 
states of neutrinos. The underlying mechanism of neu- 
trino mixing resulting in neutrino oscillations, well es- 
tablished by several experiments over the last decade 3l- 
0], is governed by the 3x3 unitary PMNS matrix §|, 
which can be parametrized using three mixing angles and 
a CP-violating phase. Evolution of neutrino flavor eigen- 
states in vacuum depends additionally on the ratio of the 
distance traveled to the neutrino energy {L/E) and the 
splitting between the squared masses of neutrino mass 
eigenstates i and j, Am^^ = m| — mf. For three neutrino 
mass eigenstates there are two independent mass split- 
tings. MINOS, a long-baseline experiment with L/E = 
0(500 km/GeV), is sensitive to the larger (atmospheric) 
mass splitting through the disappearance of muon neu- 
trinos [3]. 

The MINOS experiment uses two detectors separated 
by a distance of 734 km, both placed in the intense NuMI 
neutrino beam from Fermilab [l^l- The Near Detector 
is used primarily to characterize the NuMI beam near 
its production. The Far Detector measures the event 
rate and energy spectra after the neutrinos have trav- 
eled through the Earth's crust. In an earlier publica- 
tion MINOS presented the most precise measurement 
to date of the atmospheric mass splitting using data from 
a beam exposure of 3.36 x 10^° protons on target (POT). 
The results in this letter are based on an exposure of 
7.25 X 10^" POT, involve additional event categories, and 
employ an improved analysis methodology. 

The NuMI beam operates with 120GeV/c protons 
directed onto a graphite target of two interaction lengths. 
Positively charged liadrons produced in the target are fo- 
cused towards the beam axis by two magnetic horns. The 
neutrino beam is the product of pion, kaon, and muon de- 
cays occurring downstream of the target, primarily along 
a 675 m long decay pipe, evacuated for the first half of 
the data set but later filled with 0.9 atm helium for 
structural reasons. Data taken at different relative horn- 
target longitudinal positions and horn currents were used 
to tune the neutrino beam simulation The effect of 
the helium in the decay pipe and an observed decrease in 
neutrino flux per POT, attributed to target degradation, 
are incorporated into the simulations. Most of the data 
were taken in three run periods (Runs I— III), with the 
target placed in the most downstream position, yielding 
an energy spectrum of neutrino interactions peaking at 



3 GeV. A small amount of the data was taken with the 
target placed upstream in the high energy (HE) conflg- 
uration which results in an energy spectrum that peaks 
at 9 GeV. 

Both MINOS detectors [l3,[l3 are placed on the NuMI 
beam line axis. Each is a tracking, sampling calorimeter, 
built of 2.54 cm thick iron plates interleaved with scintil- 
lator planes composed of 1 cm thick, 4.1 cm wide scintilla- 
tor strips, arranged in two alternating orthogonal views 
and read out using multi-anode photomultipliers. The 
Near Detector is located 1.04 km downstream from the 
target and has a 23.7 1 fiducial mass. The Far Detector 
has a 4.2 kt fiducial (5.4 kt total) mass. Both detectors 
are magnetized, with a toroidal magnetic field oriented 
to focus negatively charged particles. 

In comparison to the previous analysis 0] , the data set 
has been substantially increased and several changes to 
the simulation, reconstruction, and analysis methodology 
have been introduced. This analysis also benefits from 
the inclusion of new event samples: events with a recon- 
structed track of positive charge and events originating 
outside of the fiducial volume, including the surround- 
ing rock. The geometrical modeling of the target and 
focusing system IJ] of the NuMI beam line, using the 
FLUGG software package [IlLcombining GEANT4 [l^ 
geometry with the FLUKA [17| hadron production, has 
improved the beam simulation. As previously, the detec- 
tor simulation uses GEANT3 with NEUGEN3 [Hi as the 
neutrino interaction generator. 

The most significant analysis improvement is achieved 
by employing a k-nearest-neighbor (kNN) algorithm [l^ 
for estimation of the energy of showers produced by low 
energy hadronic cascades. In MINOS, a charged cur- 
rent interaction of a muon neutrino is typically charac- 
terized by a muon track and a hadronic cascade, recon- 
structed as one or more showers. The energy resolution 
of events with a hadronic cascade is usually limited by 
the calorimetric measurement of the showers, which has 
a resolution of 56%/^E(GeV) © 2% [13 . This can be 
compared to a resolution of 4.6% or 11% for 3 GeV muon 
tracks measured by range or curvature, respectively. The 
new hadronic energy estimator uses three event charac- 
teristics: the summed reconstructed energy deposited by 
showers within 1 m of the track vertex; the sum of the en- 
ergy in the largest two reconstructed showers associated 
with the event; and the length of the longest shower [20j . 
The hadronic energy is taken as the mean true hadronic 
energy of the closest Monte Carlo events in this three 
dimensional feature space. Monte Carlo studies show 
that the new algorithm improves shower energy resolu- 
tion from 55% to 43% for showers between 1.0 — 1.5 GeV. 

The identification of charged current neutrino interac- 
tions uses the energy deposition and its fluctuation along 
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FIG. 1: The energy spectra of fuUy reconstructed events in 
the Near Detector classified as charged current interactions. 
The sohd and open circles show the data reconstructed with 
negative or positive track charge, respectively. The solid lines 
show the tuned Monte Carlo with a shaded error band due 
to systematic uncertainties. The shaded area at the bottom 
represents the simulated antineutrino component. 



a track to discriminate muon tracks from spurious tracks 
reconstructed from hadronic activity in neutral current 
interactions. However, this identification method, used in 
the prior publication 0, [2l| , does not resolve events with 
muons crossing 10 or fewer detector planes (correspond- 
ing to about 500 MeV of muon energy) . An additional 
procedure is applied to reclaim some of these events by 
constructing a new kNN classifier from the total pulse 
height of the last 5 scintillator strips associated with a 
track, and from two quantities that are correlated with 
track scattering (2^. The new selection achieves a 90% 
charged current efficiency. Below 2 GeV, the efficiency is 
77% with 6.5% neutral current contamination. 

Events classified as charged current interactions are 
used irrespective of the reconstructed charge-sign of the 
track. Compared to the previous analysis [j] , which used 
only events with a well identified negative track charge, 
events at low energy, where track charge-sign reconstruc- 
tion is less reliable, have been now recovered. Below 
6 GeV, the main oscillation energy range, 98.0% of all 
selected events arise from neutrino interactions; the an- 
tineutrino component, shown in Fig.[T] is small and con- 
tributes primarily at higher energies. This antineutrino 
component is assumed to oscillate with the same param- 
eters as the neutrinos. 

The predicted energy spectrum in the Far Detector is 
calculated from the spectrum measured in the Near De- 
tector, using a technique that takes account of the kine- 
matics of neutrino production in the beam and of the ge- 
ometry of the NuMI beam line [llj. The Near Detector 
events with tracks of positive and negative reconstructed 
charges are used separately to provide energy spectra 
predictions at the Far Detector [gl, [13]. The Far De- 



tector events with a reconstructed negative track charge 
are further divided into five quantiles based on energy 
resolution |2J| determined by simulations and test beam 
measurements [isj . This division increases the sensitivity 
with which MINOS can measure the neutrino mass split- 
ting and mixing [2^ , since events with the most precisely 
reconstructed energy carry the most precise information 
about the energy dependence of charged current event 
disappearance. 

This analysis includes interactions originating in the 
rock and outside of the Far Detector fiducial volume. 
Such interactions are only partially reconstructed, and 
are characterized b y th e measured muon and its detector 
entry position 2^ 2^. The predicted energy spectrum 
for these events is derived using the same method as for 
the fully reconstructed events. The partially and fully re- 
constructed samples have comparable statistics, but the 
partially reconstructed events contribute primarily to es- 
tablishing the overall event rate since they are due to 
neutrinos that are not well measured and are predomi- 
nantly at higher energies. 

The effect of systematic uncertainties on the measured 
oscillation parameters was determined using Monte Carlo 
simulations in which modeling parameters were varied. 
Table U shows the systematic effects, their 1 a variation 
level, and the impact on the values of mass splitting and 
mixing angle. Uncertainties in the physics simulations, 
including pion absorption cross-sections in the nucleus 
and associated modeling of energy deposition in the de- 
tector, result in the uncertainty in the visible hadronic 
energy (a), which is energy dependent and is about 7.0% 
below 3 GeV. The errors in the measurement of muon 
energy (b) from range (2% error) or from curvature in 
the magnetic field (3%) are included. The effects of rel- 
ative reconstruction efficiencies between the two detec- 
tors and uncertainties in their fiducial masses and rel- 
ative difference in detector structure result in the 1.6% 
normalization error (c). These three uncertainties dom- 
inate the systematic error on the neutrino mass split- 



Source of 


5(Am^) 


(5(sin^(2e)) 


systematic uncertainty 


(10-3 eV^ 




(a) Hadronic energy 


0.051 


< 0.001 


(b) jj, energy (range 2%, curv. 3%) 


0.047 


0.001 


(c) Relative normalization (1.6%) 


0.042 


< 0.001 


(d) NC contamination (20%) 


0.005 


0.009 


(e) Relative hadronic energy (2.2%) 


0.006 


0.004 


(f) a„(E, < 10 GeV) 


0.020 


0.007 


(g) Beam flux 


0.011 


0.001 


(h) Neutrino-antineutrino separation 


0.002 


0.002 


(i) Partially reconstructed events 


0.004 


0.003 


Total systematic uncertainty 


0.085 


0.013 


Expected statistical uncertainty 


0.124 


0.060 



TABLE I: Sources of systematic uncertainties, their one stan- 
dard deviation variation level, and their impact on fitting os- 
cillation parameters. 
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FIG. 2: Top: The energy spectra of fully reconstructed events 
in the Far Detector classified as charged current interactions. 
The dashed histogram represents the spectrum predicted from 
measurements in the Near Detector assuming no oscillations, 
while the solid histogram reflects the best fit of the oscillation 
hypothesis. The shaded area shows the predicted neutral cur- 
rent background. Bottom: The points with error bars are the 
background-subtracted ratios of data to the no-oscillation hy- 
pothesis. Lines show the best fits for: oscillations, decay [30l |. 
and decoherence [sj. 



ting. The largest uncertainty in the mixing angle is from 
the amount of neutral current background (d), the un- 
certainty on which, based on a data-driven method, is 
20% [2^. Other sources of uncertainty include: the 2.2% 
relative energy calibration uncertainty between the two 
detectors (e) ; uncertainties in the neutrino cross-sections 
Gu (f); the beam flux (g); and uncertainties due to mis- 
classification of neutrino and antineutrino interactions 
(h). Finally, incorporation of partially reconstructed 
events introduces a small uncertainty due to approxima- 
tions made in modeling the rock composition and details 
of the Far Detector's edges (i). 

All event selection criteria and analysis procedures 
were defined prior to examining the full data set in the 
Far Detector. The energy spectra were compared with 
those used in the previous publication 01 . These agree 
within the small differences expected due to changes in 
the reconstruction algorithm. The observed and pre- 
dicted numbers of events, classified in the Far Detector as 
fully and partially reconstructed charged current interac- 
tions, for all running periods are shown in Table HIl The 
energy spectrum of the fully reconstructed Far Detector 
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FIG. 3: The muon energy spectra of partially reconstructed 
events in the Far Detector. Conventions as in Fig. (2] 

data sample is shown in Fig. [2l along with the predicted 
spectra. The corresponding spectra for the partially re- 
constructed events are shown in Fig. [31 

To test the neutrino oscillation model against the data, 
the two-parameter survival probability formula P(y^ 
Va) = 1 — sin was used. The best 

values of |Am^| and sin^(2^) were found by maximizing 
a likelihood, which includes the four dominant system- 
atic uncertainties a-d in Table |I] as nuisance fit param- 
eters [H, . The likelihood value is computed at each 
point in the |Am^| — sin^(26') plane by summing the con- 
tributions from the seven event categories. Within each 
category the likelihood value is calculated by compar- 
ing the observed energy spectrum with that predicted 
for the oscillation parameters of that point. The best 
fit value and one-dimensional uncertainties for the mass 
splitting are jAm^] = (2.32l° J^) x IQ-^eV^. For the 
mixing angle, if sin^(26') is constrained to be < 1, the 
best fit is sin2(26') = 1.00 or sin2(26') > 0.94 (0.90) at 
68 (90)% confidence level (C.L.). The best fit values with 
the resulting 68% and 90% C.L. contours are shown in 
Fig. m Imposition of the physical boundary on sin^(20) 
results in the total uncertainty being smaller than that 
expected from the simulated sensitivity. Without this 



Run 


POT 


Predicted 


Observed 


Period 


(10^°) 


(No oscillations) 


(Far Detector) 






Fully 


Partially 


Fully 


Partially 


I 


1.269 


426 


375 


318 


357 


II 


1.943 


639 


565 


511 


555 


III 


3.881 


1,252 


1,130 


1,037 


977 


HE 


0.153 


134 


136 


120 


128 


Total 


7.246 


2,451 


2,206 


1,986 


2,017 



TABLE II: Numbers of events classified in the Far Detector as 
fully and partially reconstructed charged current interactions 
shown for all running periods. The predicted numbers are 
calculated under the assumption of no oscillations. 
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Laboratory, and the personnel of Fermilab for their con- 
tributions to this effort. 
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FIG. 4: Likeliliood contours of 68% and 90% C.L. around the 
best fit values for the mass splitting and mixing angle. Also 
shown are contours from previous measurements 



constraint the best fit value for mass splitting changes 
by (5(Am^) = — O.OI x IO~'^eV^ and for the mixing angle 
changes by 6ism\2d)) = +0.001. The fits do not signif- 
icantly pull away from their nominal values any of the 
four nuisance parameters. Predicted energy spectra for 
the best fit arc shown in Fig. [2] If the fit is restricted 
to use only fully reconstructed events with the negative 
track charge, the best fit value for mass splitting changes 
by (5(Am^) = -1-0.03 x 10~^eV^ and the mixing angle 
is unchanged. Two other hypotheses for neutrino dis- 
appearance, pur e neutrino decay [s^ and pure quantum 
decoherence [31|, are excluded at 7 and 9 standard devi- 
ations, respectively, as shown in Fig. [2l 

In summary, MINOS data from a beam exposure of 
7.25 X 10^° POT, more than double the data set used 
in the previous MINOS publication and improved 
analysis methodology have resulted in the measurement 
of the value of the atmospheric mass splitting to be 
Am^l =(2.32l:°og) x fO^^eV^ and the mixing angle to 
be sin2(26') > 0.90 (90%C.L.). This is the most precise 
measurement of this mass splitting to date. Neither the 
pure quantum decoherence nor neutrino decay hypothe- 
ses fit the observed spectra. 
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STFC; the US NSF; the State and University of Min- 
nesota; the University of Athens, Greece; and Brazil's 
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Minnesota DNR, the crew of the Soudan Underground 
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